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For the simple model of cold electrons streaming against cold ions the complete 
set of nonlinear stationary waves is expressed in terms of elliptic functions. The 
conditions for their dynamical connection to a uniform neutral plasma state are 
taken into account, and the conditions for the neglect of the magnetic field are 
analysed. The range of existence of stationary waves is found to be confined to 
the stable regime of the two-stream instability, but covers only part of it. All 
nonlinear BGK waves that are found within the limits of the model can be 
shown to bifurcate from the two-stream instability, some of them also 
exhibiting secondary and further bifurcations. As an exceptional case, all 
bifurcations can be treated exactly. Close to the linear regime, all nonlinear 
modes turn out to be unstable. The corresponding instability is caused by a 
wave decay that transports energy from low to high wavenumbers of the 
Fourier modes constituting the wave. From the two-stream solutions four- 
stream solutions with exactly vanishing magnetic field are derived. 


1. Introduction 


The two-stream instability and the existence of nonlinear stationary 
electrostatic waves have been well known phenomena for a long time. They 
were studied using the simple plasma model of cold electrons streaming against 
cold ions: the two-stream instability by Pierce (1948) and Bohm & Gross (1949) 
and BGK waves by Akhiezer & Lyubarskizs (1951) and Dawson (1959). More 
sophisticated plasma models were also used: a waterbag model was studied (see 
Davidson 1960 and references therein); and Bohm & Gross (1949) employed 
distribution functions for treating the two-stream instability and a simple 
example of BGK waves, while with the same method Bernstein, Greene & 
Kruskal (1957) developed a full theory of the waves named after them. A recent 
review of nonlinear electrostatic BGK modes is that of Schamel (1986). In this 
paper we return to the simple cold-plasma model and closed some gaps in the 
existing theory. The simplicity of this model is at least partly compensated by 
the possibility to obtain rather far-reaching analytical results. Of course, 
important physical phenomena like particle trapping and Landau damping are 
lost this way, but, since the BGK waves obtained for non-singular distribution 
functions are undamped, this seems to be no great drawback. Limitations that 
can be attributed to the simplifications involved in the model employed in this 
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paper will be pointed out near the end of §3.1 and in §3.5. Emphasis is put on 
the problem : where does the magnetic field vanish or where can it be neglected ? 
In analogy to solutions investigated by Holloway & Dorning (1991) within the 
framework of a plasma model with distributed particle velocities, we also derive 
four-stream solutions with two cold-ion and two cold-electron streams for which 
the magnetic field vanishes in the whole parameter regime. 

The paper is organized as follows. In §2 we discuss the basic equations of the 
problem and give a brief review of the linear theory of the basic states. In §3 
we derive a compact representation of BGK waves in terms of elliptic 
functionst{ that is very suitable for numerical purposes. An important feature 
of our approach is that we take into account all conditions arising from the 
requirement that the nonlinear waves must be dynamically connected to the 
basic two-stream state. To our knowledge, these complicating conditions have 
never so far been taken completely into account. Although the BGK solutions 
can be considered to be solitary waves on a periodic interval, it will be shown 
in §3.4 that infinitely extended solitary waves do not exist. We shall 
furthermore investigate singular solutions that represent a lattice (§3.3). In §4 
we study the connection between nonlinear BGK waves and the two-stream 
instability via bifurcations, in continuation of some previous work on this 
subject by one of the present authors (Rebhan 1983). Some preliminary results 
on the nonlinear stability of BGK waves are presented in §5. Finally, in §6 we 
extend our solutions to four-stream solutions with exactly vanishing magnetic 


field. 


2. Basic equations and linear stability of the basic two-stream flow 


Throughout this paper we shall employ dimensionless quantities. Since their 
relation to dimensional quantities will be needed for estimates, we shall briefly 
describe their relation to dimensional quantities (indicated by a tilde and 
measured in SI units). All lengths are normalized with the length 2, = 1 m, the 
densities of ions and electrons with Ža = N/&, where N is a typical number of 
electrons per cubic metre chosen such that the normalized electron density in 
the basic flow (see below) is equal to 1. The electrical field is normalized with the 
Coulomb field 
_ Ne 
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Ë (1) 
where ĉ is the electric charge of an electron and ê is the vacuum dielectric 
constant. All velocities are normalized with the velocity defined by the equation 

m_i = ah, žo, (2) 
where %%_ is the electron mass; the time is normalized with i, = #,/d, and 
magnetic fields with By = Ey/%. 

We consider perturbations of a basic flow that consists of ions with constant 
(normalized) density N, = 1 streaming with constant (normalized) velocity V, 
and of electrons with constant density N_ = 1 streaming with constant velocity 
V_in the x direction. In the perturbed states the normalized densities are given 

t After completion of this paper, one became aware of a different representation in terms 


of elliptic functions by Leven & Miiller (1974), described in an internal university 
publication. 
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by 1+n, and the normalized velocities by V,+v,. We assume that the 
problem has an intrinsic periodicity in the normalized co-ordinate x and 
decompose the normalized electric field into its period average E, which gives 
rise to a voltage drop over one period and a contribution that derives from a 


periodic potential: u 
E=E-0,®. (3) 


With the normalization described above, the flow equations for ions and 
electrons and Poisson’s equation for the coupling electric field are 


ny +0,[(V, +v,)(1+n,)] = 0, (4) 
ô v, +30,(V; #4, 18 =(, (5) 
n,—n_+0,,0 = 0, (6) 


where T = m_/Mm,. We shall assume that either the magnetic field generated by 
the plasma current and the dielectric current is exactly zero, in which case the 
equation u 
d,‚E+7=0 (7) 
(derived in the appendix) must be satisfied, or that the magnetic field must be 
so small that it can be neglected. The condition that must be satisfied in the 
latter case is 
m 10? max |n, —n_| 
d,E 0:5 —- ——_+—., 8 
EI N max|V, +04, (8) 
where X is the typical number of electrons per cubic metre (again this is derived 
in the appendix). If (7) is not satisfied but (8) is then (4)-(6) provide the so- 
called electrostatic approximation. 
In terms of the notation adopted in (4)-(6), the time-independent basic state 
is given by 
ny, =v, =0=0. (9) 
Since the equations for n, and v, are of conservation form, the mean values 7, 
and J, are conserved quantities (the averages extend over one periodicity 
length). Since all perturbed states evolve from the basic state, we must satisfy 
the conditions 
RA,=0 and 7, =0 (10) 
of dynamical connection. 
The quantities 


1 
w= gp (ttn) (V +o)? tall +n) (Ve +o)? +E”, 


(11) 
J= zn) (Veto HER (+o +d E+ 
obey the equation 0,wt+0,J = Eld, E +7), (12) 


and can be considered as the energy density and energy flux density 
respectively. @ is also a conserved quantity, provided that either Æ = 0 or 
d, E +3 = 0. 
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Ficure 1. Locus of real solutions of the linear wave dispersion relation (15) in the (0, &) 
plane. 


The dispersion relation for linear waves proportional to e“**-* has the well- 
known form 
1 T 

BR (Ber) = a 

It depends on « only through xV, and xV_, and we can shift this dependence into 

the definition of V, and V_ by setting « = 1. With this, the periodicity length of 

linear waves in 27, which will be chosen as the periodicity length of nonlinear 

waves as well. After the transformation 


õ = w—4(V, +V), R =V,- V), (14) 


1 r 


(13) assumes the form Gre Ga = 


(15) 


Four branches of real solutions exist, shown in figure 1 in an (@, Ñ) diagram. All 
four roots of (15) are real for |x| > |R|, when the two points on the branches 
where d@/dk = œ are denoted by 

R=+k&, Ö=tö,. (16) 


For —k, <<, only two roots @ are real; the other two are complex 
conjugates, one of them describing the two-stream instability. The latter sets on 
at the parameter values of the points (16). These describe marginal modes and 
mark bifurcation points as well. 

From (15) and the condition dõ/d = œ, one easily obtains 


R= HLF Gy = G(R) = A1 + PLP). (17) 
According to (14), in the original quantities we obtain marginality for 
Wy = $0, +V, +V) (18) 
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FiGURE 2. Locus of stationary (I) (---) and marginal (II) (——) modes for k = 1 in the 
(V,, V) plane (@ indicates V°’), 


and Vi-V_ = +28. {19) 


All real-w solutions proportional to e“*-* are stationary in a co-moving frame 
of reference; hence w=0 can be achieved by a Galilei transformation. 
According to (13), stationary linear modes are obtained for 


1 2 


pte i. (20) 


Stationary marginal modes are obtained by setting wa = 0 in (18), i.e. for 
= = tõ) V= FV? = F(R tö). (21) 


In the (V,, V) diagram of figure 2 they are obtained for the parameter values 
of the points where the curve (20) for stationary modes denoted by I and shown 
dashed touches the straight lines (19) for marginal modes denoted by II and 
shown full. (Note that the w, values calculated from (18) at intersection points 
of the straight lines II with other branches of curve I are non-zero, i.e. they 
belong to a different solution branch and do not describe stationary marginal 
modes.) 
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3. Nonlinear stationary-mode solutions 
3.1. One-humped regular periodic waves 
Averaging (5) with (3) over the periodicity interval yields 


T)-e 
0,04, = { je + 0. 
Therefore the equations for nonlinear stationary modes are obtained from 


(4)-(6) by dropping the time derivatives and setting Ẹ = 0. With properly 
chosen integration constants c, and W,, they can be integrated to yield 


(V,+04)(1+n,) = + (2D), | 


1 (22) 
(Vi.t+v_)(L4+n_) = +2%?, 
(V,+0,)? +210 = 2rW,, | es 
(V.+0_)2?-20 =2W., | 
rr = Nn N. (24) 


The consequences of the additional equation j = 0 or the condition (8) for the 
neglect of the magnetic field will be discussed at the end of this subsection. In 
(22) two possible signs of the integration constants have been introduced, 
admitting positive and negative signs of V}. The signs that must be chosen in 
the different quadrants of the (V,, V_) plane can be obtained from the linear case 
(v, +0,n, 0) by a continuity argument, and are shown in table 1, the first sign 
pertaining to the ions and the second to the electrons. Equations (23) can be 
solved to yield 


v, = ETW, -O)E-V,, v- = +12(W +07. (25) 


The signs to be chosen here can be concluded from the conditions (10), and 
again correspond to Table 1, where the first sign pertains once again to the ions 
and the second to the electrons. 

Inserting (25) into (22) yields 


È c? 


th. SS ee 26 
maa aE Sar (26) 
in all four quadrants, and (24) can be rewritten as 
2 2 
©" (x) +—+ _-_—_ = 0, (27) 
(W.—oF (W+o} 


Equations (25) and (26) allow the calculation of n, and v, from ®(x) once the 
latter has been determined from (27). 

Since (x) is periodic, it must have a minimum in the periodicity interval. 
The x origin can be chosen such that the minimum lies at x =0, whence 
®’(0) = 0, and, furthermore, ® can be gauged according to ®(0) = 0. With this, 
(27) can be integrated once to yield 


[D (x) ]? = 402 {(W + (a)! WE} + 4c2{(W, — (x) Wi. (28) 


Setting ©’ = 0 then yields an equation for ® that has two zeros, i.e. periodic 
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+,+ for V, >0,1” >0 —,+ for V,<0,V 
—,— for V,.<0,h <0 +,- for V, >0,M 


TABLE 1 


solutions ®(x) have just one minimum and one maximum value of ®. The 
maximum value ®,, is determined by the second zero of ®’(x) and satisfies 


(W.+0, WETO oP 
where we define Q by Q= (=). (30) 
From (2), we get ®(x) = +2c, W(®(x)), (31) 
where W(z) = {QUW +2)! Wi] + (W, —2)t— WEË, (32) 


i.e. ®’(x) has two branches. ®(x) is symmetric around its extrema, and for the 
one-humped 2rr-periodic solutions considered in this subsection the maximum 
is attained at x = 7. 

Equation (31) can be integrated once more to yield the implicit solution 


(2) dz 


i We = 2c, x (33) 


for the first branch of (x) extending from the minimum to the maximum, and 
a corresponding one for the second extending from the maximum to the 
minimum. Solutions ®(x) whose first branch would extend from maximum to 
minimum (®’(x) < 0) can be obtained from the latter ones as ©(x) = ®(z+) 
owing to the translational invariance of (27). 

It is possible to express the implicit solution (33) for ® in terms of elliptic 
functions. For this purpose, we first transform from W, to new parameters 


W: 
a,=—, (34) 
+ Ou 
in terms of which Q becomes 
4 a _4\ 
ga ee (35) 
(a_+1)!—a? 
We now go from z = ® to a new variable y by setting 
_ siny 
siny = TA (36) 
where u = arcsin (=u) a, u, = arcos f (37) 
ta 
3 1 
with igs Pee ee. (38) 
Q [(a, +a) (1 +Q°)} 
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The lower limit of integration z, = 0 in (33) transforms to u, = arccos f and 
yw, =4n. Since sinu < sinu, for z > 0 according to (37), y is defined for all 
z < ®,,. From (32) and (36)-(38), we get 


W(z) = Aib yla, +a_)(1+Q*)} sin du, cosy (39) 
and 
dz MOy(a,+ayFf R  _ 2 sin? YH — ksin? yy) 
We a+ Facer, 8Qk? sin? y(1— k? sin? y) 
+210°-1)(ksiny- sin’ |dy, (40) 
where k is defined by k = sintu.. 


Since cosu, = 1-2sin’iu, = 1—2k? = 2, we get, with (38), 


1 at +Q 


k=l — Atv l 41 
al aon et 


y in {m y 
Using f | ay...=-| ay’...+ f dy’ ..., 
Yı y 0 0 
fe sin? y'(1—k?sin? y’ idy’ = 41 — k?) f eae A 
0 o (1—ksin? y} 


+4(2k?— 1) i (1— k? sin? ydy’ —}k? sin y cos y(1— k? sin? y}, 
0 


and the notation 
E(k, p) = | : (1—ksin? y’)idy’, F(k, yw) = [= (42) 
0 o (1— k? sin? py’) 
for the elliptic functions, from (33) with (40) and with the abbreviations 
E=EK(k,in), F =F(k,in), 
we finally get the implicit representation of our solution ®(z): 
(2k? — 1) E(k, yr) — (k? —8) F(k, wr) — k? sin y cosy(1— k? sin? y)? 


3k(Q?— 1) 


_ 92 4 2h one? 
+ FTA) (1—2k? + 2k? cos? y) cos y 


= (2k —1)E-(k ante” (43) 


Then, according to (36)-(38), ® = z is related to y by 


Pula, ta) 


= 
1+@ 


{1 — 2k? + 2kQ(1 — hk)? 


—[1—2k? sin? y+ 2kQ sin (1 —k? sin? yt}. (44) 
With the solutions thus obtained, we have to satisfy the four conditions (10) 
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and the periodicity condition ® = ®,, for x = x, = 7. According to (36)-(38), x, 
transforms tö Y, = —}7, and, inserting this in (43), with E(k, —}7) = —E and 
F(k, —im) = —F, we obtain 


3nc,(1+Q?) 


k? —3) F — (2k? —1)E = +. 45 
Ppr ee-e ak (45) 
According to (26) with (31), the condition 7, = 0 becomes 
$ =I 
ens 2n| meal l (46) 
0 W(z2)(W,— 2) 
With 
dz MO la, +a_)} QU —k sin? y)—2ksin y(1— k? sin? y) dy 
W) W, =z} AdE} (1—k? sin? y}? 
this can be rewritten as 
2 3 
n(1+Q¥ (47) 


+ = 9Q2H —F)[O,(a, to) 


The condition 2_ = 0 is automatically satisfied if 7, = 0, because averaging (24) 
over one periodicity interval yields 2_ = 7,. Eliminating c, from (45) and (47) 
yields 

377(1 +Q?)? 


Pu = ra, +a.) OEP (eF Oe —1) By 


(48) 


We have finally to satisfy the conditions (10), which with (25) can be rewritten 


as 
_ (T$ [Pm (W,—2)} _ 2 [0 (W_+z} 
m= + m [ 2c, W(z) oe =| 2c, Wey” er 


Since V, and V_ do not enter the solution ®(x) when parametrized by a, the last 
two conditions can be neglected for the moment. 

Discussing the existence of solutions in the (&,,&_) plane now becomes very 
simple. From the five conditions to be satisfied by our solutions, only the first 
three relate ®,, and the integration constant c, to the parameters a, and a_. 
These parameters enter our solutions via the quantities Q and k’, and the latter 
exist and are real if a! and («,—1)! are real, i.e. for «_ > 0 and a, > 1. The 
physical meaning of these conditions is the exclusion of particle trapping, 
because in the interior of the regime of existence they imply 


(V,+4,)? = 2T(W,-®) > 2T@,,(a,—-1) > 0 
and (V_+v_)* = 2(W_+®) > 20,,a_ > 0 for all x. The boundaries of existence in 
the (a,,@_) plane are given by the two straight lines 
a_=0, a,=1. (50) 


In summary, solutions can be obtained in the following way. We start by 
prescribing values a_ 20 and a, > 1. For these, we calculate Q from (35), k? 
from (41) and c, from (47). This way, all parameters are given that are needed 
to determine our nonlinear solution ®(x) from (43) and (44). The amplitude ®,, 
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FIGURE 3. ®,, as a function of a_ and a,. 


is shown as a function ofa_ and a, in figure 3. The largest amplitude is obtained 
for a_ = 0 and a, = 1, and is given by 


_ In? 
SHE? (1+2) EF +41+2)F?] 


where E and F must be evaluated for k? = }(1 —274), Once ®(x) is determined, 

we get n, from (26) after we have calculated W from (34) and c? from (30). 

Equations (49) can then be used to determine the original parameters V, that 

correspond to a,. Once these are determined, we finally obtain v, from (25). 
V, and V_ can also be expressed in terms of elliptic functions. Using 


Ox x 4638 ..., 


(W,—2)i = el [1 — 2k? sin yy + 2kQ sin y{l—k?sin? yh, 
(W. +2) = | [Q(1 — 2k? sin? y)—2ksin y(1— k? sin? WH], 


F k? sin? y (1— k? sin? ydy = (1 — k?) F + (2k?— 1) E], 
tn 


f ksin? y (1— ksin? ydy” = 4[(1 — k?) (2k?— 1) F + (4kt—2k?—1) E], 
-gr 


and (40), after some calculations we obtain 


2r „la, + = Ps BQO, (x, +a_)]* 


v,=+ : 
2 +| 1+Q? 15rc,(1+Q?} 


(0+ Gp) B+ ar + Qbr) FI}, 


E T e] l _ Aody(a,+a_)] (51) 


2 2 
1+Q* 1570, (1 +02} (art+br)E+(Q apt bp) FI}, 
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ap =—96kt+56k?—1, bp = 32k*—32k'4+2, 


h 
Mises ap = 48kt—34k?+1, bp = —16k9+ 18k?—2. 


Our result (51) can be used to transform the region of existence for nonlinear 
solutions into the original (V,, V_) plane. Figure 4(a) shows the existence regime 
that was obtained this way in the fourth quadrant of the (V,,V_) plane. 
According to (68) and table 1, the existence regime in the other quadrants is 
obtained from this by simple reflections at the V, and V_ axes respectively. The 
boundary corresponding to a. = 0 (or a, = 1 respectively) was calculated by 
inserting pairs a_ = 0,a, > 1 (or a_ 20,a, = 1 respectively) and the value ®,, 
following from them in (51). An additional bounding curve is provided by the 
linear-mode limit ®,, 0. It is obviously represented by (20), and corresponds 
to infinity in the («,,a_) plane according to (34). Basically, the curves of 
existence for linear solutions are broadened into band-like regions of existence 
for nonlinear solutions that are confined to the stable regime of the two-stream 
instability. Apparently the latter restriction is closely related to the loss of 
particle trapping in our plasma model. 

The curves used for shading the region of existence are curves with constant 
a_> Q or constant a, > 1 respectively. It is seen that they cover some parts of 
it more than once. Thus, in contrast to the situation of the (a,,a_) plane, we 
have regions with multiple solutions in the (V,, V_) plane. Figure 4(b) presents 
a schematic view of these regions for the fourth quadrant, and indicates the 
multiplicity of the solutions. 

In figure 5 we show what the solutions ®(x) and the corresponding densities 
1-+n,(x) look like at various points of the existence regime marked 1-6 in figure 
4 (a). We see that the densities approach a sinusoidal shape (point 1) at the left- 
hand boundary of the existence regime, which is due to the limit of linearity, 
and develop more and more pronounced peaks (points 2-6) towards the right- 
hand boundary. The latter is due to the appearance of singularities. According 
to (34), a. = 0 and/or a, = 1 imply W_ = 0 and/or ®,, = W, respectively, and 
it follows from (26) and (27) that our solutions become singular then. More 
precisely, n_(x) and/or n,(x) become infinite (but remain integrable), whereas 
(x) remains regular. Depending on the point of the right-hand boundary that 
is approached, the profile steepening can occur for electrons (point 3) and ions 
(point 4) separately or for both simultaneously (point 6). 

We come now to the consequences of the so-far neglected equation for the 
magnetic field. Let us first consider the case where the latter vanishes exactly. 
In this case we must have 


j= (L+n,)(V,+0,)—(L+n_)(V+0_) = 0, 


which reduces to V, = V_ for the limiting case of linear waves. With (25) and 
(26), the condition becomes 


Q= 
in the first and third quadrants of the (V,, V_) plane. In accordance with the case 
of linear waves, the condition Q = —I® obtained for the second and fourth 


quadrants cannot be satisfied, since, according to its definition (30), Q must be 
positive. Figure 6 shows the regions of existence of nonlinear solutions in the 
first and fourth quadrants of the (V,, V_) plane and the curve Q = T? where the 
magnetic field vanishes exactly. Furthermore, for given values of the typical 
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Figure 4. (a) Existence regime of nonlinear solutions in the fourth quadrant of the (V,, V) 
plane. (Similar existence regimes can be obtained for the other three quadrants by simple 
reflections.) The curves used for shading the existence regime are curves with constant a_ > 0 


Bifurcation of BGK waves 13 


number N of electrons per cubic metre, it shows the regions where the magnetic 
field can be neglected according to (8). Corresponding domains in the second 
and third quadrants are obtained from these by reflection at the origin. 
Allowing for an error of 3x 10°? (i.e. B<3x10°E,,,„/max(V,+v,)), with 
j = Z(+ VT-Q) and the choice $,,/7 instead of max |n, —n_| for the typical 
electric field, the boundaries of these regions are obtained from the equations 


2 r(Q+T?) (52) 


max |V, +v] 10" 
— ttt - +4. 
Oy N 
(On the left-hand side, for the first/second and third/fourth quadrants of the 
(V,,V_) plane the —/+ sign must be taken.) In the second and fourth 
quadrants, owing to Q < 0 and the + sign on the left-hand side, (52) can only 
be satisfied for the + sign on the right-hand side, and the boundary is only 
obtained for N < 101!/1-77. 

In §6 we shall construct four-stream solutions from our nonlinear two-stream 
solutions for which the magnetic field vanishes exactly in the whole regime of 
existence. 

3.2. Multi-humped regular periodic waves 

Solutions with n maxima and minima in the periodicity interval [0, 27[ can be 
obtained from the one-humped solutions discussed so far by a simple 
transformation. We take advantage of the fact that we have shifted the 
dependence of our solutions into the velocities V, by setting « = 1. Replacing 
V, > V,,/mand rescaling the functions n ,, v, and ® properly, we obtain solutions 
with mean velocities V} /m and m maxima in the original periodicity interval of 
length 27. If n(x), v} (x) and ®(x) with mean velocities V, are solutions of 
(4)-(6) then the rescaled functions 


(Kx), Kku (kz), K°O(Kx) (53) 

corresponding to the rescaled parameters 

KV,, KW, Ric, 
are also solutions of (4)-(6) with periodicity length 27x. Choosing x = 1/m, we 
get a solution with m maxima and minima in the interval 27. One- and multi- 
humped waves with equal mean velocities are only obtained when both V, and 
V,,/m lie in the existence regime depicted in figure 4. This will not be the case 
in general. 

3.3. Lattice solutions 
In addition to the solutions obtained in §3.2, which are regular or limiting cases 
of regular ones, there exist isolated singular solutions for V, = V_ = 0 that will 


be considered now. Equations (22)-(24) together with the side conditions (10) 
can also be solved using Dirac delta distributions. It is easily seen that the 


function a (0<x<n), 
wi eo (7 <x < 2n), oe) 
n_=2n6(x)—1, n, =2n6(a—7m)—-1, (55) 
+ (21) '[©,,—®(x)}—const, v_ = +[20(x)}!—const, (56) 


or constant a, > 1. For the points marked 1-6 the profiles O(x), 1+m_(x) and 1 +n,(x) are 
shown in figure 5. (b) Schematic view of the regions in (a) with different solution 
multiplicities, the latter being indicated by the numbers in each region. 
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2n 
Figure 5. Solutions ®(x) (——), 1+ n_(x) (------ ) and 1+n,(x) (---) for the points numbered 
1-6 in the existence regime of figure 4(a): (1) V, = 0:05, V. = —1-16; (2) 0-04, — 1-45; (3) 0-05, 
— 1-64; (4) 0-037, — 2:2; (5) 0-041, 7:8; (6) 0-046, — 1-96. 


0 es Li. _ on 
0-6 ( 


)-8 1-0 1-2 1-4 


Fiaure 6. Domains of negligible magnetic field in the (V,,V_) plane: the magnetic field 
vanishes exactly on the dotted line, while in the shaded regions 
B < 3x 10°F „/max (V,+v,). 


In the first quadrant (52) yields two boundaries of the shaded region, those shown 
corresponding to N= 10" electrons m”®. In the fourth quadrant only one boundary is 
obtained, that shown corresponding to N = 3:3 x 10™ electrons m=’. 


is a solution in the distributional sense. From the maximum value of 
D(x), Py = n°, it may be seen that this solution is not connected to the previous 
ones. It corresponds to the case where the electrons and ions are at rest in delta- 
distribution peaks that are regularly spaced on a lattice. Of course, each peak 
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represents a two-dimensional plane that extends perpendicular to the lattice 
direction. It is easily verified that, in the same way, one can construct further 
lattice solutions with more and irregularly distributed peaks in one periodicity 
interval. 


3.4. Non-existence of infinitely extended solitary waves 


The existence of singular waves in a periodic interval would lead to regular 
solitary solutions of infinite extension if in (53) the spread of the interval were 
combined with an adequate reduction of the wave amplitude. This, however, is 
not the case. In this subsection we show quite generally that infinitely extended 
waves with a finite number of humps do not exist. 

There are two possibilities for solutions of this kind: 


0-0, as |z| >œ or @®>0 as |z|> o. 


In both cases ©’(0) +0 and ®’(0)—0 as |x|> 00 are necessary conditions, the 
second being satisfied automatically according to (31) and (32). From (27), (30) 
and (34) for ®’(0) = 0, we get 


_ (%_+0/Oy i 
ee) ne 
In the first case (PD > ®,,) this leads with (35) to the condition 
a —(a.,—1)' _ (art) 
(a_+1)i—at (a,—1)} 


(58) 


which is satisfied only on the line a_+a, = 0 in the (a_,a,) plane. In the second 
case (P>0) we get from (35) the condition 


al — (a, 1) _ (=). (59) 
(a_+1)F—at \a, 

which leads again to a_+a, = 0. The states satisfying this condition lie outside 
the range of existence (a_ > 0,a, > 1) for stationary waves. 

The non-existence of infinitely extended solitary waves must be considered as 
a drawback of the plasma model employed in this paper, because such waves do 
exist in plasmas with distributed velocities (Schamel 1986). Since they always 
come along with trapped particles, the non-existence in our model is probably 
due to the lack of the latter. 


3.5. The direction of wave propagation 
The equation V, = (%,V,+m_V)/(m,+m_) = (V+ TV )/(1+T) for the 
plasma velocity in dimensionless quantities transforms to the same equation 
Vp = (V,+TV_)/(1+1) in dimensionless quantities. Since V, and V_ are the 
velocities of ions and electrons respectively in the rest frame of the waves, Vp 
is the plasma velocity in this frame. Correspondingly, 


HABA 


Jy = "TAT 


(60) 


is the velocity of the waves in the rest frame of the plasma. Waves with zero 
velocity are obtained on the straight line V, = —I'V_, which, owing to the 
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smallness of T, almost coincides with the V, axis of the (V,, V_) plane. To the 
right of this line all waves run to the left, and to the left of it all waves run to 
the right. The existence of nonlinear waves running in both directions provides 
the possibility for studying the collision of nonlinear waves. 


4. Bifurcation of stationary nonlinear solutions 


We shall now show that the nonlinear states derived in the last section 
bifurcate from the basic states (9) at the curve (20) that represents stationary 
linear modes. (For brevity, we shall restrict our analysis to the case of one- 
humped waves.) This way, we obtain one of the rare examples where for a 
rather complex system the bifurcation of states can be followed exactly over the 
full range of existence. Although we can use our nonlinear results directly for 
the treatment of the bifurcation problem, we shall nevertheless also briefly 
sketch the steps of a bifurcation analysis, because it allows us to follow the 
solutions directly in the space of the original parameters (V,, V_). 

A genuine bifurcation analysis has to start at the critical point given by (21), 
and searches for nonlinear solutions on some curve that leads through this 
point. A possible choice for this curve would be a straight line, and the direction 
of it could be chosen such that it runs through any point selected in the (V_, V,) 
plane. For the moment, we shall consider only the straight line 

V.=V84+7, W=V-r (61) 
parametrized by 7 that runs under 45° from the left down to the right in the 
(V,, V_) plane. (It corresponds to the horizontal line © = @, in the diagram of 
figure 1.) 

Since the nonlinear solutions ®(x) can be chosen to be symmetrical with 
respect to x = 7, a Fourier expansion of ®(x) with cosine terms only is possible. 
It turns out that a bifurcation expansion of the form 

D(x) = ecos x + ep” cos 27+ 62O® cos 3x4... (62) 
with the small amplitude factor e is successful when inserted into (22)-(24). 
(Note that, with this ansatz, ®’(0) = 0 automatically, so that the expansion 
agrees with our treatment of the full nonlinear case in the previous section.) It 
is seen from (25) that n, (x) and v, (x) have the same symmetry as ®(x), whence 
again an expansion using only cosine terms is possible, and an ansatz of the 
form 


f(a) = ef cosa+ ef cos 2x + ef) cosatf% cos 3x)+... (63) 
turns out to be successful for all of them. Obviously, with this ansatz, the side 
conditions (10) are satisfied automatically. Finally, the integration constants 
contained in (22)-(25) and the driving parameter 7 in (61) that describes the 
distance from the critical point must also be expanded in terms of e, the ansatz 

Bench W= Wet... TERE... (64) 
being successful this time. Inserting all the series expansions into (22)-(24), to 
first order in e we get the flow equations 


n@V? +P =0, 


pry © =0, 
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with solutions 


T 1 
M D= a a a 
no 7= o? nl a? Vy + yo’ vt nE ye (65) 
and the electric equation n®—n®+09,,00=0, (66) 
which, with (62), leads to E i 1 (67) 
3 > yo ye 


as a solvability condition. Comparison with the condition (20) for stationary 
linear solutions shows that a bifurcation in the direction of the straight line (61) 
is possible not only from the critical point (21) but also from any point on curve 
I in figure 2. 

To order e?, we obtain the equations 


r 
(2) y0 (2) L Jap (1) „1 — (2) V0 4 1,,(1)? (2) _ 
nen Hol Hi =0, oP VE HN +{ "ho = 0, 


nD — nD —40® = 0, 
with solutions 


(68) 
et T 3m mil 31 
* AVE 4ye 48 ap 
(vr 1 
p? = mw) (69) 


The direction of the bifurcation along the straight line (61) is determined by the 
sign of r,, which can only be obtained at third order in e. There, the flow 
equations yield solutions that depend on 7,, and inserting them in the third- 
order electric equation yields an equation for 7, with solution 


_ 1 (TV F2/(V9)t + (11/V2 —3)/(V2)* 


= 3 3 70 
2776 r/ve—1/V® en 


Since V_ is negative, 7, is positive for all points on the curve for linear stationary 
modes. Thus the bifurcation leads everywhere into the region where the two- 
stream instability is stable. According to (64), we have a pitchfork bifurcation, 
which means that for any wave amplitude e, —e is also possible, i.e. with 
P(e, x), P{-—e,x) is also a solution. The second solutions ®(—e,x) = ®(e,c-+ m) 
are obtained from the first ones by a trivial translation. Figure 7 shows a 
bifurcation diagram for the bifurcation along the straight line (61). The main- 
mode amplitude e obtained from the bifurcation analysis is shown by the 
upper clear straight line, while the one extracted by Fourier analysis of the full 
nonlinear mode is shown by the line below this, with dots at the calculational 
points. The bifurcation turns around nowhere, and ends at the boundary of 
existence W_ = 0. It is thus confirmed that in the unstable regime of the two- 
stream instability no stationary nonlinear modes exist. Obviously each point on 
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Figure 7. The upper two curves show the amplitude of the first cosine mode along the line 
(61), obtained from bifurcation analysis ( ) and from the exact nonlinear solution 
(--— — ). The lower two curves show the difference between the total energy of the 
bifurcating states and the corresponding basic state obtained from bifurcation analysis 
( ) and from the exact nonlinear solution —+—+—). 


©, 35 
3-0 
2-5 
2-0 
1-5 
1-0 
05 


Ficure 8. Secondary bifurcation along a straight line in the (V,, V_) plane crossing the 
region where three nonlinear solutions exist (cf. figure 4b). 


the boundaries of the region in the (V,, V_) plane with two or three nonlinear 
solutions (see figure 4b) is a possible critical point for a secondary bifurcation. 
For example, if we give the bifurcation line (61) a different inclination such that 
it runs into the region with three solutions, we get a first bifurcation at its point 
(V?, V°?) defined by (21) and a second one just where it crosses the boundary of 
this region (figure 8). 


5. Stability of the nonlinear solutions 


In the last section we saw that the nonlinear solutions bifurcate from the 
basic state into the stable region of the two-stream instability. Under certain 
conditions, one can conclude from this behaviour that in the vicinity of the 
critical point the bifurcating states are unstable and that the system tries to 
return from them to the corresponding basic state (same values of V, as the 
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bifurcating states, but n, = 0 and v, = 0). This conclusion is not possible here: 
numerical calculation shows that the total energy W = f = w(x) dx contained in 
one periodicity interval of the bifurcating solutions is always smaller than the 
total energy of the corresponding basic state (see figure 7). The same result was 
obtained from the bifurcation analysis, yielding 


= = ani Vase. V_—2c, E} —2c_Ei = W,,,+.a(V9, V°) e4+O(e*), (71) 
where the energy of the basic flow is Was = V2/2T +4V? and where a(V°, V?) 
stands for a very lengthy expression that has negative values for the bifurcating 
states. From this, we can conclude that the system cannot try to lower its 
energy — provided that some friction is present — by returning to the basic state, 
and thus it is not possible to draw any conclusions from the bifurcational 
behaviour regarding stability. 

So far we have carried out only a very simple stability analysis that does not 
allow very far-reaching conclusions. The quantities ®(x), n,(x) and v, (x) were 
expanded in Fourier series 


fee) 
F(x) = 3 (a cosvata sinve), F=®,n,,v,, (72) 


v=0 


etc. and inserted into (4)-(6). The coefficients of ® were eliminated, yielding a 
finite system 


Ne), w=s,e, v=0,...,N, 


of coupled nonlinear amplitude equations after truncation at some finite N. This 
system was integrated numerically with a fourth-order Runge-Kutta method. 
The start values of the amplitudes were chosen close to an exact nonlinear 
solution, i.e. at start time t = 0 the amplitudes a% were equal to zero and the 
a were calculated by Fourier expansion of a numerically obtained nonlinear 
solution. Perturbations of the exact state came about by numerical uncertainty 
and/or a small uniform reduction or enhancement of all Fourier coefficients. 
Many nonlinear solutions were tested, and in the regime of small nonlinearity 
all turned out to be unstable. The reason for this instability was found to be a 
transfer of wave energy to large mode numbers. This seems to indicate either 
turbulent wave decay or a transition into a time-dependent state with strongly 
peaked waves, and will be investigated in greater detail in a forthcoming study. 


6. Periodic four-stream waves with exactly vanishing magnetic field 


In this section, from the two-stream solutions considered so far, we derive 
stationary nonlinear four-stream solutions with two ion and two electron 
streams for which the magnetic field vanishes exactly in the whole parameter 
regime. For this purpose, we assume one stream of ions/electrons with density 
N, +a, at velocity V,+v,, and for each species a second one with density 
N,; Hna, at velocity —(V, +v,), such that the basic flow satisfies the condition 


N,,+N,,=1, +N, =1. (73) 


20 L. Hannibal, E. Rebhan and C. Kielhorn 
The equations corresponding to (22)-(24) are 

(Vi. +v) (Mi +744) = r), (74) 
—(V, +05) (No, tn) = Ere, (75) 
(V, +v,)’+2T0 = 2TW,, (76) 
(V_+v_)?—20 = 2W,, (77) 
(78) 
where r, =T, T_=1. (79) 


Ong PB = Ni- Hng — N14 Nor 


The combination of signs must again be chosen according to table 1. The 
magnetic field vanishes exactly for 


J = (Nis +14 — Nog — nog) (Va +04) — Ni- +i- Na- —n,_)(V_ +o) 
= Ae}, — 03) +(e —e8_)] = 0, (80) 


where the + sign holds in the second and fourth quadrants and the — signs in 
the first and third quadrants of the (V,, V_) plane. 
From (74)-(77), we obtain 


cÈ, ce 
= — N ; mi = —*__-N, 3 = 1,2, 81 
"+ won oo "aoe eh (81) 


and inserting this in (78) yields 


e, +c c_+c2 
te re, 82 
mol (+o) a 


For g +t, =, ete =e (83) 


this reduces to (27), and we obtain from (80) the additional condition 


2 _ 22 ; 
a (84) 
1-7 - 


(with + in quadrants 1 and 3 and — in quadrants 2 and 4 of the (V,, V_) plane). 
From any two-stream solution of (27), we obtain a current-free four-stream 
solution by choosing ¢,,, ¢,,, ¢,_ and c,_ such that the three nonlinear equations 
(83) and (84) are satisfied. This is always possible, because, with the ansatz 


Cz =CzSINPy, Cy = Cy COSP,, 


(83) are satisfied and, using (30), (84) becomes 


cos2¢, _ +2 
cos2d_ “TE 


Since, by proper choice of ¢, and ¢_, the left-hand side can assume all possible 
values between — œ and œ, this equation can always be satisfied. 


This research was performed under the auspices of the Sonderforsch- 
ungsbereich 191 of Deutsche Forschungsgemeinschaft. 
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Appendix. Conditions for vanishing or negligible magnetic field 


If the magnetic field vanishes exactly, the Maxwell equation that relates the 
magnetic field to the currents becomes 


0,£+j=0, where j= (1+mn,)(V.t0,)—(L+n_)(L+v_), (85) 
in our dimensionless quantities. Since from (3), (4) and (6) we get 


‘x m x 
=f m-ni- f aa | dx’ (n’,—n_)+E, 
27 Jo ü 


0 
with the consequence 6,£+j =d, E +3, 
(85) is satisfied for all x if 


d,‚E+7=0. (86) 
For the stationary solutions considered in this paper this condition reduces to 
j=0. (87) 


When (87) is not satisfied, the current j will generate a magnetic field B. In that 
case we consider a plasma cylinder with radius r whose axis extends in the x 
direction and we assume cylindrical symmetry of the magnetic field. Instead of 
(86), we then get 


~\2 
B=ta(d,H+j)r, with a= (=) , (88) 


where B is the poloidal magnetic field at radius r and £ is the velocity of light. 
The magnetic field would have to be taken into account in the momentum 
equations (5). It is strictly negligible if 


IV, +04) Bl < |El, (89) 


for all x, and it can be neglected on the average if the inequality holds for a 
typical Ep, which can be obtained from (6) as 


Eryp 
With (88), r x 1 and a = 35N x 107" 


typ? 


x max |n,—n_|. 


according to (1) and (2), where N is the typical number of charges per cubic 
metre, the condition (89) for the neglect of the magnetic field finally becomes 


10™ max |n, —n_| 


ld, E +71 < ET aCe 


(90) 
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